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ABSTRACT 

Context. The helioseismic determination of the solar age has been a subject of several studies because it provides us with an inde¬ 
pendent estimation of the age of the solar system. 

Aims. We present the Bayesian estimates of the helioseismic age of the Sun, which are determined by means of calibrated solar 
models that employ different equations of state and nuclear reaction rates. 

Methods. We use 17 frequency separation ratios ro 2 (n) = (v„,/=o - Vn-i,i=i) I - v„_i,; = i) from 8640 days of low-? BiSON frequen¬ 

cies and consider three likelihood functions that depend on the handling of the errors of these r 02 (n) ratios. Moreover, we employ the 
2010 CODATA recommended values for Newtons constant, solar mass, and radius to calibrate a large grid of solar models spanning a 
conceivable range of solar ages. 

Results. It is shown that the most constrained posterior distribution of the solar age for models employing Irwin EOS with NACRE 
reaction rates leads to f E = 4.587 ± 0.007 Gyr, while models employing the Irwin EOS and Adelberger, et al.. Reviews of Modern 
Physics, 83, 195 (2011) reaction rate have t Q = 4.569 ± 0.006 Gyr. Implementing OPAL EOS in the solar models results in reduced 
evidence ratios (Bayes factors) and leads to an age that is not consistent with the meteoritic dating of the solar system. 

Conclusions. An estimate of the solar age that relies on an helioseismic age indicator such as ro 2 (n) turns out to be essentially 
independent of the type of likelihood function. However, with respect to model selection, abandoning any information concerning 
the errors of the ro 2 (n) ratios leads to inconclusive results, and this stresses the importance of evaluating the trustworthiness of error 
estimates. 

Key words. Sun: helioseismology - Sun: interior - Sun: oscillations - Equation of state - Nuclear reactions, nucleosynthesis, abun¬ 
dances - Methods: statistical 


1. Introduction 

By definition, the age of the Sun is the time since the proto¬ 
sun arrived at some reference state in the HR diagram, usually 
called thejjbirth line” where it begins its quasi-static contraction 
(Stable). 1 19831) . For a star of 1 Af Q , this stage depends on var¬ 
ious factors, such as rotation, spin, and wind accretion during 
the protostar phase, and it should produce a pre-main-sequence 
(PMS) object o f ~ 4000 K with a luminosity of ~ 10 L 0 
dstahler & Pallall200l . 

The age of the solar system is instead assumed to lie between 
the age of crystallized and melted material in the solar system, 
and the time of significant injection of nucleosynthesis material 
in the protosolar nebula. The latter is an upper limit to the solar 
age, and the precise relation between the “age of the Sun” and 
that of the planetary bodies and nebular ingredients depends on 
the detailed physical conditions during the collapse of the proto¬ 
solar nebula. 


Recent studies based on the dating of calcium-aluminum- 
rich inclusions (CAI) in chondrites have reported a mete¬ 
oritic “zero age” of the solar system r a nging roughly from 


4.563 to 4.576 Gyr 

Bahcall et al., 19951c Amelin et all 120021 

Jacobsen et all 2008, 

2009; Bouvier & Wadhwa, 20lll). A more 


recent estimate supports an age of 4.567 Gyr ( Connelly et al. . 
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12012 ). The interesting question is to locate this time in the PMS 
evolution of the Sun as defined 

To what extent tis he solar system age consistent with the no¬ 
tion of “birth line” and its location in the HR diagram? In fact, 
the PMS evolution is often neglected in stellar evolution calcu¬ 
lation, based on the duration of this phase being much shorter 
than the successive evolution. On the other hand, the precise lo¬ 
cation of the zer o-age ma in-se quence (ZAMS) is problematic, as 
discussed in lMorel et al. : ( 2000 ). 

Several studies have thus tried to use helioseismology to pro¬ 
vide an independent estimation of the solar age, thus testing the 
consistency of the radioactive dating of the solar system with 
stellar evolution theory. The standard approach to calibrating he¬ 
lioseismic solar age has been to confronting specific oscillations 
diagnostic among solar models of different ages, while keeping 
the luminosity, th e radius, and the heavy elements abundance at 
the surface fixed ( Gough & Novotnvl 19901: Dziembo wski et all 


1999: Bonanno et al. , 20021 Christensen-Dalsgaardl 2009t) . The 

resulting helioseismic age is then a “best-fit” age that can 
be obtained with a series of calibrations. The limitations 
of the one -param eter calibration approach ( Gou sfl HqoI 
iHoudek & Goughl 1201 lh he in the difficulties of estimating the 
effect of ch emic a l compositi on and, in general for unknown 
physics dBonanno et all 200l|), of determining the sound speed 
gradient near the center (lGoughl 12012 ). 

The essential ingredient for estimating the helioseismic so¬ 
lar age is to find a “genuine” oscillation diagnostic for which 
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the sensitivity to the complex physics of the oute r layers of 
the star is minimal. As i s well known dRoxburgh & Vorontsovi 
12003b lOtf Floranes et all 120051) . the frequency separation ratio 
n,i+ 2 (n) = (V„,/-V„-u+ 2 )/(>«,/+ 1 —V,,-!,/+ 1 ) of low-degree p-modes 
is a fairly good indicator of the inner structure of the sun because 
it is mostly sensitive to the gradient of mean molecular weight 
near the center. In particular in lDo&an et al.l 1 201(f ). it was shown 
that ro 2 (n) is a relatively robust age indicator since it does not de¬ 
pend on thesurface-effect correction of the higher order p-modes 
dKieldsen et all 120081) or on different definitions of the solar ra¬ 
dius. 

In this work we discuss a Bayesian approach to the he¬ 
lioseismic determination of the solar age. In fact in recent 
times, the use of the Bayesian inference has bec ome increas¬ 
ingly common in the astronomical community ( Trottal 120081) . 
In the case of asteroseismology, the use of Bayesian infer¬ 
ence has been essential for establishing credible and robust 


ina (Ouirion et all 2010: Bazot et al. 

, 20l2 

Gruberbauer et aQ 

20fj). In Gruberbauer & Guenther 

2013b 

t has been argued 


that that there is an inconsistency between the meteoritic and the 
solar age as inferred from helioseismology. In our case the ad¬ 
vantage of using the Bayesian inference as opposed to the stan¬ 
dard frequentist’s approach lies in the possibility of rigorously 
comparing the relative plausibility of solar models with differ¬ 
ent physical inputs, chemical compositions, etc., by computing 
evidence ratios (i.e., Bayes factors). 

Besides comparing models with different EOS and nuclear 
reaction rates, we use the new recommended 2010 CODATA val¬ 
ues for the solar mass, the Newton constant, and the solar radius, 
and we check the robustness of our findings by using different 
likelihood functions to quantify the agreem ent of model predic¬ 
tions with the BiSON data (iBroomhall et al.lf2009i) . We show that 
updated physical inputs lead to an helioseismic age that is con¬ 
sistent with the meteoritic age of the sol ar system, a s also sug¬ 
gested from the analysis of lGruberbauer & Guenther] d201 3l) . On 
the other hand, at the level of accuracy of current helioseismol¬ 
ogy, it is essential to include the PMS evolution (about 40-50 
Myr) in the standard solar model calibration. 

The structure of this work is the following. Section 2 com¬ 
piles the physical details of the solar models considered and the 
observations. Sect. 3 describes the Bayesian approach, Sect. 4 
presents the results, and Sect. 5 is devoted to the conclusions. 


2. Models physics and observations 

Our non-rotating solar model s were built with t he Ca tania ver¬ 
sion of the GARSTEC code (iWeiss & Schlattl [20081) . a fully- 
implicit ID code that includes heavy-elements dif fusio n. It em¬ 
ploys either the OPAL 2005 equation of state (iRogers et al.LI 1 99(1 
iRogers & NavfonovU2002i) . compl emented with the MHD equa - 
tion of state at low temperatures (iHummer & Mihalasl Il988l) . 
or the Irwin equation of state (ICassisi et all 2003b, and it uses 
OPAL opacities for high temperatures (ilglesias & Rogersl 1996 ) 
and Ferguson’s opacities for low temperatures (Ferguson et al, 
120051) . The nuclear reaction rate s are either fro m the NACRE col¬ 
labora tion (lAngulo et all 1 19991) or from the lAdelberger et alJ 
(1201 ll) co mpilation, and the chemic al composition follows the 
mixture of iGrevesse & Noelsl (i 19931) with (Z/X) Q = 0.0245 at 
the surface. We also consider models with th e so -called “new 
abundances’’ for which (Z/X) Q = 0.0178 ( Asplund et a l., 120091 ). 

Our starting models are chemically homogeneous PMS 
models with log L/L B - 0.21 and log T e = 3.638 K, so they are 
fairly close to the birth line of a 1 M Q object, which according 


to IStahler & Pallal d2005i) . would be located about 2 Myr 
before. The value of Newton’s constant is the 2010 CODATA 
recommended value G = 6.67384 x 10~ 8 cm 3 g _1 s _2 taken from 
http://physics.nist.gov/cuu/Constants/index.html 
As a co nsequ ence, because GM 0 = k — 1.32712440 x 10 26 
cm 3 s” 2 ( Cox , 2000) , we assumed M 0 = 1.98855 x 10 33 g. The 
radius is taken to be R e = 6.95613 x 10 10 cm based on an 
average of the two valu es and quot ed error bar in Tab le 3 of 
iHabeireiter et al.l (i2008l) . (See also iReese et aD (1201 2l) for an 
application of the 2010 CODATA to seismic inversions.) The solar 
luminosity is instead L e = 3.846 x 10 33 erg s -i (ICoxll2000l) . 

We then used the definitive “best possible estimate” of 
8640 days of low-/ frequency BiSON data, corrected for 
the solar cycle modulation dBroomhall et all 120091 ) taken from 
http://bison.ph.bham.ac.uk/index.php?page=bison,frequei 
In particular, we considered N - 17 frequency separation ratios 
roiiri) = (v n ,/ =0 - v n _i,/ =2 )/(v n ,/=i - v„_ u= i), ranging from order 
n — 9 to order n — 25 for the / = 0,1,2 modes, together with the 
corresponding uncertainties. 


3. Bayesian inference 

In the Bayesian view the central quantity to be computed is a 
conditional probability P(A\B): It represents the probability that 
event A will happen given the fact event B has actually occurred. 
In our case it is a subjective type of probability, a degree of plau¬ 
sibility of A given B, and as such it bears no resemblance to a 
frequency distribution. In common applications A is a vector of 
parameters that quantify a model and B the set of observational 
data. Estimates of the parameters have the same conceptual sta¬ 
tus as probabilistic events. 

According to Bayes’ theorem, 

P(A\B) = P(A)^^ (1) 

holds, where P(A) and P(B) are unconditional probabilities for 
events A and B. In our application P(A\B) is the probability that 
model parameter r falls within the interval r... r + dr given the 
data. P(B\A) is the usual likelihood function, P(A) the so-called 
prior for r, and P(B) the searched-for evidence of the model. 
The evidence is a model’s mean likelihood - averaged over the 
whole r range and subject to the prior probability density - and 
as such, it measures the overall ability of the model to cope with 
the data. Since the prior probability must sum up to unity, the 
evidence diminishes if the parameter space expands beyond the 
space needed to cover the essential parts of the likelihood moun¬ 
tain. 

We are interested here in evidence ratio, the Bayes factors. 
They allow a model ranking, i.e. to compare the explanatory 
powers of competitive models or hypotheses. The logarithm of 
age, r = log e (f), is chosen, with t being the solar age, to make 
certain that the posterior for the age is compatible with the pos¬ 
terior of, say, the reciprocal of the age. With this decision the ev¬ 
idence does not depend on the unit of age. Accordingly, for P(A) 
a flat prior is assumed over the logarithm of age. All the other 
parameters of the solar model - the initial helium fraction To, 
the mixing length parameter, and the initial (Z/X) o value - are 
adjusted in order to reproduce measured radius, effective temper¬ 
ature, and Z/X ratio at the surface. We do not attempt to derive 
credibility intervals for To, (Z/X)o, and the mixing length pa¬ 
rameter, because they are obtained analytically, i.e., by means of 
the standard Newton-Raphson procedure embedded in the solar 
model calibration. 
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In actual calculations the likelihood function P(B\A ) is often 
expressed by a Gaussian. In particular, if d, = r^Jn) are the 
observed data (n = i + 8, i = 1.. .N), m, the theoretical model 
values, and cr,- the errors, it reads as 


N 1 

ET 1 

(d, - m,(r)) 2 * 

r— eX P 
i=l V27rcr,- 

2cr 2 


( 2 ) 


with r being log e (f) and N = 17. An observed quantity, d, = 
ro 2 (/+ 8), is the ratio of two Gaussians that results, strictly speak¬ 
ing, in a Cauchy-like distributiorQ not a Gaussian. With respect 
to the following likelihood function, which is based on Eq. 0, 
and for the sake of consistency, we decided, however, to com¬ 
pute the cr,-’s according to Gauss’s linearized error propagation 
rule. 

Although [Broomhall et al.l ( 2009 ) have carefully taken the 
systematic error induced by solar activity into account, it is 
tempting to enhance the robustness of the age estimates by av¬ 
eraging over the errors, thereby obeying Jeffreys’ 1 /cr prior. In 
this way an error-integrated likelihood is obtained, which relies 
still on the Gaussianity assumption: 


a' = r°A— = 

Jo cr 


T (N/2) 


2 (nN) N/2 


[n; v =, -v,| i2f =1 (^)' 


N/2 


(3) 


where cr, = ,v, • cr, and relative errors ,y, are normalized according 
to Wi = N, with Wi - 1/s 2 the weights. 

We also consider the median likelihood, which is constructed 
from the binomial distribution and ignores any information with 
regard to measurements errors. It only assumes that both positive 
and negative deviations dj - m, have equal probability. It reads 
as 


A" 


n 

;=i 


1 N\ 

2v Kl(N-K)'. ’ 


(4) 


with K the number of events with d, > m, (or d, < mi). 

The evidence is the mean likelihood over parameter space; 
that is to say, in the case of a flat prior for r = log e (f). 


1 


Tu - T1 



A (,!) dr, 


(5) 


where the interval t u - ri is the same in all cases where, when 
considered otherwise, the communicated Bayes factors would 
render useless. 


4. Results 

To properly resolve A(r), the r domain was covered by 
80 equidistant grid points. We checked that our results 
are substantially insensitive to a further refinement of the 
grid. (The achieved log time resolution is, if expressed in a 
musician’s language, a little bit better than one cent.) The 
computation of the eigenfrequencies has been performed 
with the latest version of the public pulsation code, GYRE 

1 One can indeed substitute Eq. m by the correct expression. The 

corresponding ages and their standard deviations are indistinguishable 

from the values communicated in Table [7] However, as the distribution 

is a “fat-tailed” one, all Bayes factors but one (that normalized to unity) 

are somewhat enhanced: in the worst case (“Irwin+AdelR+Asplund”) 

by 44 per cent, in all other cases by up to 14 per cent. More important, 

the ranking is not affected. 


Table 1. Expectation value and standard deviation of the solar 
age for the Gaussian likelihood function (0 with known errors 
cr,. The ordering of the models is according to their Bayes fac¬ 
tors. 


Input physics 

UGyr] 

Bayes factor 

Irwin +NACRE 

4.587 ± 0.007 

1 

Irwin + AdelR 

4.569 ± 0.006 

0.17 

OPAL +NACRE +old CODATA 

4.696 ± 0.006 

2.1 x 10~ 13 

OPAL +NACRE 

4.702 ± 0.006 

1.9 x 10^ 15 

OPAL + AdelR 

4.683 ± 0.008 

3.0 x 10- 17 

Irwin +AdelR+Asplund 

4.785 ± 0.006 

2.9 x 10~ 44 


Table 2. Ages and Bayes factors in the case of unknown magni¬ 
tude of the over-all error cr in Eq. 0. 


Input physics 

GlGyr] 

Bayes factors 

Irwin +NACRE 

4.587 ±0.010 

1 

Irwin +AdelR 

4.569 ±0.010 

0.42 

OPAL +NACRE +old CODATA 

4.697 ± 0.017 

2.3 x 10“ 4 

OPAL +NACRE 

4.703 ±0.018 

1.0 x 10~ 4 

OPAL + AdelR 

4.683 ±0.018 

5.6 x 10- 5 * 

Irwin +AdelR+Asplund 

4.786 ± 0.026 

1.1 x 10~ 7 


Table 3. Ages and Bayes factors in the case of median statistics, 
Eq. 0. The ordering of the models is as in Tables 0 and [2] 


Input physics 

fe[Gyr] 

Bayes factors 

Irwin ±NACRE 

4.582 ±0.014 

0.38 

Irwin ±AdelR 

4.563 ±0.015 

0.43 

OPAL ±NACRE +old CODATA 

4.692 ± 0.026 

0.60 

OPAL ±NACRE 

4.698 ± 0.028 

0.63 

OPAL ±AdelR 

4.676 ± 0.029 

0.65 

Irwin ±AdelR± Asplund 

4.775 ± 0.039 

1 


https://bitbucket.org/rhdtownsend/gyre/wiki/Home, 
which employs a n ew Magnus multiple-shooti ng scheme, as de¬ 
scribed in detail in iTownsend & Teitlerl d2013l) . For consistency 
we have thus implemented the 2010 CODATA also in GYRE . The 
results, ages, and Bayes factors are depicted in Tables 00 for 
the three likelihood functions 00 considered. 

First, we note that the use of a significantly longer data 
set dBroomhall et all 1200 9). in combination with a Gaussian 
likelihood 0, has provided us with a much sharper sola r ag e 
e stimat e if comp ared to previous studies (I Do van et al.L 12010 : 
IChristensen-Dalsgaardl 12009 ): models obtained with Irwin EOS 
generally perform better than model s with OPAL EOS , while the 
evidence for models with NACRE and Adel berger et al.l (120 ll) re- 
action rates differ only marginally. Only Adelberger et al.l ( 2011 1 
reaction rates, combined with an Irwin EOS result in a helioseis¬ 
mic age t Q , are consistent with the meteoritic one within lcr. 

The m ain reason for the age diff erence between models with 
NACRE and I Adelberger et akl (1201 ll) reaction rates is the value of 
the astrophysical .S’ pp (0)-factor for the pp-fusion cross-section. 
The NACRE collaboration adopts S pp = 3.89x10 25 MeVb, while 
[Adelberger et al.l (1201 ll) adopt Spp = 4.01 x 10~ 25 MeVb. As 
a consequence, the helioseismic age tends to be longer for the 
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NACRE rates because of the reduced efficiency of the pp reac¬ 
tions. 

Using the 2010 CODATA does not produce a significant im¬ 
pact on the solar age determination. On the other hand, if we 
consider for the OPAL +NACRE case in Table 0, which is the 
old CODATA value for Newton’s constant G = 6.67232 x ICE 8 
cm 3 g _1 s' 2 , we find t Q = 4.696 ± 0.006 Gyr and a Bayes factor 
that is two orders of magnitude larger. This value for the solar 
age is consistent within one lcr with the value of f Q = 4.62+0.08 
Gyr found by iDogan et al.l (2010) (including PMS evolution), 
which was obtained with a much shorter BiSON data set. 

In the case of the error-integrated likelihood (0, the expec¬ 
tation values of the solar age are reported in Table [2] They are 
basically the same as in the former (Gaussian) case (Table Q]). 
but the uncertainties are at least a factor two larger. Of course, 
this is due to the extreme stance that even the magnitude of the 
error cr in 0 is assumed to be unknown. 

Even more robust ages estimates are presented in Table [2 
In median statistics, nothing is assumed about the shape of the 
symmetric error distribution. The expectation values of the solar 
age are systematically lower by up to 10 Myr as compared with 
the the two other cases where the Gaussianity of the errors is pre¬ 
sumed. The credibility regions are up to seven times larger than 
in the case of a Gaussian likelihood with trusted errors. More 
important, in terms of Bayes factors, all models now perform 
equally well! 

5. Conclusions 

Bayes factors are powerful indicators when it comes to quan¬ 
tifying the ability of solar models differing in input physics to 
cope with published “quiet Sun’’ frequency separation ratios. If 
one trusts the common Gaussianity assumption, models using 
the Irwin EOS perform best. Abandoning any error information, 
to be on the safe side, leads to inconclusive results with respect 
to Bayes factors (cf. Table 0, which stresses the importance of 
evaluating the trustworthiness of the error estimates that enter 
into Eqs. 0and0. 

If we assume a Gaussian error distrobution, only models 
with Irwin EOS agree with the meteoritic age. Moreover, their 
Bayesian evidence exceeds those with an OPAL EOS by at least 
four orders of magnitude (cf. Table 0. With the Adelberger et 
al. reaction rates, the solar age proves to be even more consistent 
with the meteoritic age. 

Incorporation of PMS evolution is essential because the stan¬ 
dard deviation of the age estimation, < 10 Myr in the case of 
Irwin EOS, is less than the 40-50 Myr time span of the PMS 
phase. Moreover, our helioseismic age is consistent with the no¬ 
tion of a birth line, since our starting PMS models are only 2 
Myr distant from the birth line location. The switch from old 
CODATA to 2010 CODATA values does not effect the age esti¬ 
mate significantly, in contrast to the model’s evidence, at least if 
the published errors are taken seriously. 

These conclu s ions are obtained using the “old” 
iGrevesse & Noelsl (1 1993b with (Z/X ) 0 = 0.0245 at the 

surface. It is worth mentioning that the so-c alled “ne w abu n- 
dances’,’ for which ( Z/X) Q = 0.0178 ( Asnlund et al .. 2009 '). 
would lead to completely inconsistent values of the solar age 
and Bayes factors reduced by many orders of magnitude as 
indicated by the corresponding row entries in Tables [7] and [2] 


calcium-aluminum-rich inclusions and the meteoritic age of the 
solar system. 
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